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ENDOCARDIAL MAPPING SYSTEM 

CROSS-REFERENCE 
A catheter suitable for use in implementing this invention i 
further disclosed in copending U.S. Patent Application serial 
number, not available, filed, not available. 


TECHNICAL FIELD 

The present invention relates to a system for measuring 
electrical signals originating within cardiac tissue and generating 
a map of the electrical activity of the heart. 

BACKGROUND ART 

It is common to measure the electrical potentials present on 
the interior surface of the heart as a part of an 
electrophysiologic study of a patient's heart. Typically such 
measurements are used to form a two-dimensional map of the 
electrical activity originating from the heart. An 
electrophysiologist will use the map to locate centers of ectopic 
electrical activity occurring within cardiac tissue. 

One traditional mapping technique involves a sequence of 
electrical measurements taken from mobile electrodes inserted into 
the heart chamber. The electrodes are repetitively moved into 
contact with the endocardial surface. 

An alternative mapping technique takes essentially 
simultaneous measurements from a floating electrode array to 


generate a two-dimensional map of electrical potentials. This non- 
contact technique is taught by a number of references including 
Taccardi U.S. Patent 4,649,924. 

The two-dimensional maps of the electrical potentials at the 
endocardial surface generated by these traditional processes suffer 
many defects. One defect is related to "spatial averaging". 
Traditional systems have been limited in resolution by the number 
of electrodes used. The number of electrodes dictated the number 
of points for which the electrical activity of the endocardial 
surface could be mapped. Therefore, progressions in endocardial 
mapping have involved either the introduction of progressively more 
electrodes on the mapping catheter or improved flexibility for 
moving a small mapping probe with spot electrodes from place to 
place on the endocardial surface. Both "ring" shaped and "spot" 
shaped electrodes are constrained by spatial averaging within the 
blood volume in the heart chamber. The ring shaped electrcds 
spatially averages the measurement of electrical activity around 
the circumference of the electrode site. The spot electrode also 
spatially averages the measurement of electrical activity within 
the conical view of the electrode site. Direct contact with 
electrically active tissue is required by most/systems in the prior 
art in order to minimize this limitation of circumferential spatial 
averaging of the ring electrode or conical spatial averaging of the 
spot electrode. When a ring shaped electrode is in contact with 
the electrically active tissue, the local electrical potential at 
the point of contact is directly coupled to the electrode and the 


contact potential dominates over the effect of circumferential 
spatial averaging of the electric field. Thus the problem of 
spatial averaging is the inability to accurately resolve the 
location of ectopic tissue masses. In the prior art, iso- 
potentials are interpolated and plotted on a rectilinear map which 
can only crudely represent the unfolded interior surface of the 
heart. Such two-dimensional maps are generated by interpolation 
processes which "fill in" contours based upon a very limited set 
of measurements. They are not three-dimensional and thus cannot 
be used to distinguish a large ectopic center located a long 
distance from an electrode site, from a smaller ectopic centers 
located closer to the electrode site. The traditional two- 
dimensional mapping process also cannot locate infarcted tissue 
which is electrically inactive. The inability to accurately 
characterize the size and location of ectopic tissue frustrates the 
delivery of certain therapies such as "ablation". 

SUMMARY DISCLOSURE OF INVENTION 

The invention discloses the apparatus and method for forming 
a true three-dimensional electrical map of the interior the heart 
chamber, and a related technique for forming a two-dimensional 
subsurface map at a particular location in the endocardial wall. 

According to the present invention an intra-cardiac multi- 
electrode mapping catheter assembly 14 is inserted into the heart. 
The mapping catheter assembly 14 includes a multi-electrode array 
catheter 20 with an integral reference electrode or, preferentially 
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a companion reference catheter 16. The array catheter 20 is 
deployed to form a substantially spherical array 19 of electrode 
~~ sites. This electrode array 19 is spatially referenced to a point 

on the endocardial surface by the reference electrode or catheter 
5 16 which is moved into contact with the endocardial surface 18. 
3 The mapping catheter assembly 14 is coupled to interface 

apparatus 22 which contains a signal generator 32, and voltage 
acquisition apparatus 30. Preferably, in use, the signal generator 
32 is used to measure the volumetric shape of the heart chamber 
through impedance plethysmography. This signal generator is also 
used to determine the position of the reference electrode within 
the heart chamber. Next, the signals from all the electrode sites 
on the electrode array 19 are presented to the. voltage acquisition 
apparatus 30 to derive a three-dimensional, instantaneous high 
resolution map of the electrical activity of the entire heart 
chamber volume. This map is calibrated by the use of a surface 
electrode 24. The calibration is both electrical and dimensional. 
Lastly this three-dimensional map, along with the signal from an 
intramural electrode 26 preferably at the tip of the reference 
catheter 16, is used to compute a two-dimensional map of the 
intramural electrical activity , within the wall. The two- 
dimensional map is a slice of the heart wall and represents the 
subsurface electrical activity in the heart wall itself. 

Both of these "maps" can be followed over time which is 
desirable. The true three-dimensional map also avoids the problem 
of spatial averaging and generates an instantaneous, high 
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resolution map of the electrical activity of the entire volume of 
the heart chamber and the endocardial surface. This three- 
dimensional map is an order of magnitude more accurate and precise 
than previously obtained interpolation maps. The two-dimensional 
map of the intramural slice is unavailable using prior techniques. 

BRIEF DESCRIPTION OF THE DRAWINGS 
The various figures of the drawing represent in illustrative 
embodiment of the invention, but it should be understood that many 
modifications to the invention may be made without departing from 
the scope of the invention. Throughout these drawings identical 
reference numerals refer to identical structure throughout in 
which: 

FIG.l is a schematic view of the system. 

FIG. 2 is a schematic view representing the display of the 
three-dimensional map. 

FIG. 3 is a schematic view representing the display of the 

subsurface two-dimensional map. 

FIG.4A is a schematic flow chart of the steps in the method. 
FIG.4B is a schematic flow chart of the steps in the method. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Hardware Description 

FIG.l shows the mapping system 10 coupled to a patient's heart 
12. The mapping catheter assembly 14 is inserted into a heart 
chamber and the reference electrode 24 touches the endocardial 


surface 18. 

A preferred catheter is disclosed in companion copending U.S. 
Patent Application serial number, not available filed, not 
available, which is hereby incorporated by reference. However, any 
of a variety of catheter constructions can be used for carrying out 
the invention. For purposes of this disclosure the important 
characteristics of the catheter assembly 14 are as follows. 

The preferred array catheter 20 carries at least twenty-four 
individual electrode sites which are coupled to the interface 
apparatus 22. The preferred reference catheter 16 is a coaxial 
extension of the array catheter 20. This reference catheter 16 
includes a surface electrode site 24 and a subsurface electrode 
site 26 both of which are coupled to the interface apparatus 22. 
It should be understood that the electrode sites 24 and 26 can be 
located directly on the array catheter. The array catheter 2 0 may 
be expanded into a known geometric shape, preferably spherical. 
Resolution is enhanced by larger sized spherical shapes. A balloon 
or the like should be incorporated under the electrode array 19 to 
exclude blood from the interior of the electrode array 19. The 
spherical shape and exclusion of blood are not required but they 
materially reduce the complexity of the calculations required to 
generate the map displays . 

The reference electrode 24 or 26 and/or the reference catheter 
16 serves several purposes. First they stabilize and maintain the 
array 19 20 at a known distance from a reference point on the 
endocardial surface 18 for calibration of the shape and volume 


calculations. Secondly, the surface electrode 24 is used to 
calibrate the electrical activity measurements of the endocardial 
surface 18 provided by the electrode array 19. 

The interface apparatus 22 includes a switching assembly 28 
which is a multiplexor to sequentially couple the various electrode 
sites to the voltage acquisition apparatus 30, and the signal 
generator apparatus 32. These devices are under the control of a 
computer 34. The voltage acquisition apparatus 30 is preferably a 
12 bit A to D converter. A signal generator 32 is also supplied 
10 to generate low current pulses for determining the volume and shape 
of the endocardial chamber using impedance plethysmography, and for 
determining the location of the reference catheter. 

The computer 34 is preferably of the "workstation" class to 
provide sufficient processing power to operate in essentially real 
15 time. This computer operates under the control of software set 
forth in the flow charts of FIG. 4A and 4B. 

Software Description 

The illustrative me..hod may be partitioned into 2 6 steps as 
20 shown in FIG.4A and 4B. The partitioning of the step-wise sequence 
is done as an aid to explaining the invention and other equivalent 
partitioning can be readily substituted without departing from the 

scope of the invention. 

At step 41 the process begins. The illustrative process 
25 assumes that the electrode array assumes a known spherical shape 
within the heart chamber, and that there are at least twenty- four 
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electrodes on the electrode array 19. This preferred method c,r» 
be modified to accommodate unknown and non-reproducible, non- 
spherical shaped arrays. The location of each of these electrode 
sites on the array surface must be known. A method ot determining 
the location of the electrode array 19 and the location of the 
heart chamber walls (cardiac geometry) must be available. Since 
the physical movements of the heart are quite slow relative to the 
electrical activity there is no need to sample the geometry during 
every process cycle. A typical sampling rate for geometry update 
would be 60 Hz which is much less than the 500 Hz required to 
measure electrical activity. Therefore step 42 asks whether the 
proper interval has transpired for updating geometry measurements. 
If updating is required, then step 43 is performed. For the 
measurements in step 4 3 the preferred embodiment is the use of the 
technique of impedance plethysmography. The availability of the 
electrode array for plethysmography makes this method preferred. 
Another method for determining geometry is the use of intra-cardiac 
ultrasound to determine the catheter and the wall locations. I., 
either case, if the reference catheter 16 is extended to the 
chamber wall 18 then its length can be used to calibrate the 
geometry measurements since the calculated distance can be compared 
to the reference cathe-r length and the geometry calculation 
forced to converge on the known spacing represented by the physical 
catheters. In an alternative embodiment reference electrode 24 is 
positioned on array catheter 20 and therefore its position would 
be known . 


in step 44 the signals from all the electrode sites in the 
electrode array 19 are sampled by the A to D converter in the 
voltage acquisition apparatus 30. These measurements are stored 
in a digital file for later use in step 55. At this point (stP P 
45) the known locations of all the electrodes on the electrode 
array 19 and the measured potentials at each electrode are used to 
create the intermediate parameters of the three-dimensional 
electrical activity map. This step uses field theory calculations 
presented in greater detail below. The components which are 
created in this step are stored in a digital file for later 

use in steps 48, 52, 57 and 51. 

At the next decision point (step 46) the question is asked 
whether the reference catheter 16 is in a calibrating position, 
in the calibrating position, the reference catheter 16 projects 
directly out of the array catheter 20 establish a length from 
the electrode array 19 which is a known distance from the wall 18 
of the heart chamber. This calibration position may be confirmed 
using fluoroscopy. If the catheter is not in position then the 
process moves to step 51. This step is not required with reference 
electrode 24 positioned on array catheter 20. 

If the reference catheter 16 is in the calibrating position 
then in step 4 7 the exact position of the reference catheter 16 is 
determined using the distance and orientation data from step 43. 
In step 48, the available information includes position in space 
of the reference catheter 16 on the chamber wall 18 and the 
intermediate electrical activity map parameters of the three- 
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dimensional map. Using these two sets of information the expected 
electrical activity at the reference catheter surface electrode 
site 24 is determined. In step 49 the actual potential at this 
site 24 is measured from the reference catheter by the A to D 
converter in the voltage acquisition apparatus 30. Finally in step 
50 a scale factor is adjusted which modifies the map calculations 
to achieve calibrated results. This adjustment factor is used in 
all subsequent calculations of electrical activity. 

At step 51 the system polls the user to display a three- 
dimensional map. If such a map is desired then in step 52 a method 
of displaying the electrical activity is first determined. Second 
an area, or volume is defined for which the electrical activity is 
to be viewed. Third a level of resolution is defined for this view 
of the electrical activity. Finally the electrical activity at all 
of the points defined by the display option, volume and resolution 
are computed using the field theory calculations and the adjustment 
factor mentioned above. These calculated values are then used to 
display the data on computer 34 . 

FIG. 2 is a representative display 71 of the output of process 
57. In the preferred presentation the heart is displayed as a wire 
grid 36. The iso-potential map is overlaid on the wire grid 36 and 
several iso-potential lines such as iso-potential line 38 are shown 
on the drawing. Typically the color of the wire grid 36 and the 
iso-potential lines will be different to aid interpretation. The 
potentials may preferably be presented by a continuously filled 
color-scale rather than iso-potential lines. The tightly closed iso- 


potential line 39 may arise from an ectopic focus present at this 
location in the heart. In the representative display 71 of process 
52 the mapping catheter assembly will not be shown. 

in step 53 a subthreshold pulse is supplied to the surface 
electrode 24 of the reference catheter 16 by the signal generator 
32. in step 54 the voltages are measured at all of the electrode 
sites on the electrode array 19 by the voltage acquisition 
apparatus 30. One problem in locating the position of the 
subthreshold pulse is that other electrical activity may render it 
difficult to detect. To counteract this problem step 55 starts by 
subtracting the electrical activity which was just measured in step 
44 from the measurements in step 54. The location of the tip of 
the reference catheter 16 (i.e. subsurface electrode 24), is found 
by first performing the same field theory calculations of step 45 
on this derived electrode data. Next, three positions in space are 
defined which are orthogonally positioned near the heart chamber 
walls. The potentials at these sites are calculated using the 
three-dimensional electrical activity map. These potentials are 
then used to triangulate, and thus determine, the position of the 
subthreshold pulse at the surface electrode 24 of the reference 
catheter 16. If more accurate localization is desired then three 
more points which are much closer to the surface electrode 24 can 
be defined and the triangulation can be performed again. 

At step 56 the user of the system determines whether the 
reference catheter surface electrode 24 is in the desired location. 
If not then in step 57 this position in space can be displayed by 


superimposing it on the map of electrical activity used in step 52. 
An example of such a display 71 is presented in FIG. 2. 

W:.en step 58 is reached the surface electrode 24 is in a known 
position on the endocardial surface 18 of the heart: chamber which 
is proper for determining the electrical activity of the tissue at 
that site. If the intramural or subsurface electrode 26 which 
preferentially extends from the tip of the reference catheter 16 
is not inserted into the tissue then in step 59 the user of the 
system extends the subsurface electrode 26 into the wall 18. In 
step 60 the potentials from the surface electrode 24 and from the 
intramural subsurface 26 electrode are measured by voltage 
acquisition apparatus 30. In step 61 a line 21 along the heart 
chamber wall which has the surface electrode 24 at its center is 
defined by the user of the system. The three-dimensional map 
parameters from step 45 are then used to compute a number of points 
along this line including the site of the reference catheter 
surface electrode 24. These calculations are adjusted to conform 
to the measured value at the reference catheter surface electrode 
24. in step 62 a slice of tissue is defined and bounded by this 
line 21 (FIG. 2) and the location of the intramural subsurface 
electrode 26 (FIG. 3) and computed positions such as 23 and 25. 
Subsequently a two-dimensional map of the electrical activity of 
this slice of tissue is computed using the center of gravity 
calculations detailed below in the section on algorithm 
descriptions. Points outside of the boundary of the slice cannot 
be computed accurately. In step 63 this map of electrical activity 


within the two-dimensional slice is displayed as illustrated in 
FIG. 3. in this instance the iso-potential line 17 indicates the 
location within the wall 18 of the ectopic focus. 

In step 64 the system user's input is queried as tc whether 
the procedure is finished or not. If so it is ended (step 66). 
:f not the cycle repeats (step 65) at 500 Hz or higher, reentering 
at step 42. 

Description of the Preferred Computing Algorithms 

Two different algorithms are suitable for implementing 
different stages of the present invention. 

The algorithm used to derive the map of the electrical 
activity of the heart chamber employs electrostatic 
volume-conductor field theory to derive a high resolution map of 
the chamber volume. The second algorithm is able to estimate 
intramural electrical activity by interpolating between points on 
the endocardial surface and an intramural measurement using center 
of gravity calculations . 

In use, the preliminary process steps identify the position 
of the electrode array 19 consequently the field theory algorithm 
can be initialized with both contact and non-contact type data. 
This is one difference from the traditional prior art techniques 
which require either contact or non-contact for accurate results, 
but cannot accommodate both. This also permits the system to 
differentiate small distributions of electrical activity at the 
endocardial surface from large distributions of electrical activity 
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underlying the same surface position in the myocardium. 

in the first algorithm, from electrostatic volume-conductor 
field theory it follows that all the electrodes within the solid 
angle view of every locus of electrical activity on the endocardial 
surface are integrated together to reconstruct the electrical 
activity at any given locus throughout the entire volume and upon 
the endocardium. Thus as best shown in FIG. 2 the signals from the 
electrode array 19 on the catheter 20 produce a continuous map of 
the whole endocardium. This is another difference between the 
present method and the traditional prior art approach which use 
the electrode with the lowest potential as the indicator of cardiac 
Z]l abnormality. By using the complete information in the algorithm, 

the resolution of the map shown in FIG. 2 is improved by at least 
a factor of ten over prior methods. Other improvements include: 
the ability to find the optimal global minimum instead of 
sub-optimal local minima; the elimination of blind spots between 
electrodes; the ability to detect abnormalities caused by multiple 
foci; the ability to distinguish between a localized focus of 
electrical activity at the endocardial surface and a distributed 
path of electrical activity in the more distant myocardium; and the 
ability to detect other types of problems including ischemic 
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tissue 


The algorithm for creating the 3D map of the cardiac volume 
takes advantage of the fact that myocardial electrical activity 
instantaneously creates potential fields by electrotonic 
conduction. Since action potentials propagate several orders of 


magnitude slower than the speed of electronic conation, the 
potential field is quasi-static. Since there are no significant 
charge sources in the blood volume, Laplace's Ejuation for 
potential completely describes the potential field in the blood 
volume: 

V*<t> = 0 

From electrostatic field theory, the general sph3ricai 
harmonic series solution for potential is: 

♦ (r.8.»> - g^lV * Bl r-^}<l> la Y la (6^) (i) 

in spherical harmonic's, Y^S.V) is the spherical harmonic 
series made up of Legendre Polynomials. is the lm th component 

of potential and is defined as: 

where V(0,<p) is the measured potential over the probe radius R and 
dQ is the differential solid angle and, in spherical coordinates, 
is defined as: 

dQ = sinddBd? 

During the first step in the algorithmic determination of the 
3D map of the electrical activity each * la component is determined 
by integrating the potential at a given point with the spherical 
harmonic at that point with respect to the solid angle element 
subtended from the origin to that point. This is an important 
aspect of the 3D map; its accuracy in creating the 3D map is 
increased with increased numbers of electrodes in the array and 
with increased size of the spherical array. In practice it is 
necessary to compute the » to components with the subscript 1 set to 


4 or greater. These * u components are stored in an 1 by in array 
- : for later determination of potentials anywhere in the volume within 

the endocardial walls. 

The bracketed expression of equation 1 (in terms of A lf B lf 
5 and r) simply contains the extrapolation coefficients that weight 
U the measured probe components to obtain the potential components 

anywhere in the cavity. Once again, the weighted components are 
summed to obtain the actual potentials. Given that the potential 
is known on the probe boundary, and given that the probe boundary 
10 is non-conductive, we can determine the coefficients A, and B,, 
yielding the following final solution for potential at any point 
within the boundaries of the cavity, using a spherical probe of 

radius R: *<r.e.*) - ± £ * (^r i ]*^<e.,>) 

i-o a>--i 

15 The integral for *^ can be evaluated with an estimating sum, 

discretized by p latitudinal rows and q longitudinal columns of 
electrodes on the spherical probe. 

Note that p times q equals the total number of electrodes on the 
20 spherical probe array. The angle * ranges from zero to n radians 
and V ranges from zero to 2n radians . 

At this point the" determination of the geometry of the 
endocardial walls enters into the algorithm. The potential of each 
point on the endocardial wall can now be computed by defining them 
25 as r, 6, and <fi . During the activation sequence the graphical 
representation of the electrical activity on the endocardial 
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surface can be slowed down by 30 to 40 times to present a picture 
of the ventricular cavity within a time frame useful for human 
viewing . 

A geometric description of the heart structure is required in 
order for the algorithm to account for the inherent effect of 
3 spatial averaging within the medium. Spatial averaging is a 

function of both the conductive nature of the medium as well as 
the physical dimensions of the medium. 
3 Given the above computed three-dimensional endocardial 

potential map, the intramural activation map of FIG. 3 is estimated 
by interpolating between the accurately computed endocardial 
potentials at locations 23 and 25 (FIG. 2), and actual recorded 
endocardial value at the surface electrode 24 and an actual 
recorded intramural value at the subsurface electrode 26 site. 
This first-order estimation of the myocardial activation map 
assumes that the medium is homogenous and that the medium contains 
no charge sources. This myocardial activation estimation is 
limited by the fact that the myocardial medium is not homogeneous 
and that there are charge sources contained within the myocardial 
medium. If more than one intramural point was sampled the 
underlying map of intramural electrical activity could be improved 
by interpolating between the endocardial surface values and all the 
sample intramural values. The center-of-gravity calculations can 
be summarized by the equation: 

m - ^ 

s ^ - ^ 
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